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Abstract-This paper deals with the analysis and simulation of a 
Tunisian Wind farms in Sidi-Daoud, the most important wind 
farm in Tunisia. In this farm constant speed wind turbines are 
used. The first section devoted to the presentation of the wind 
farm of Sidi-Doud. In the second part, we present the 
aerodynamic model of the wind turbine and the transient model 
of the induction generator used in this wind park. Finally, we 
simulate the electrical behavior of the wind farm under the 
turbulent wind upstream of each rotor turbine and identify the 
main features based on simulation software PSAT. 


Keywords— Constant speed wind turbine, Transient model of 
induction machine; Power System software PSAT 


Nomenclature 


H „ : The turbine inertia constant 

H „ : The generator inertia constant 
K: Shaft stiffness 

D, : The mutual damping 

D,,, : The turbine self-damping 

Dp : The generator self-damping 

T „ : The aerodynamic torque 

T, : The electrical torque 

W : The rotational speed of the rotor 
w, : Speed of the output shaft 

w,, : The rotational speed of the generator 


sear : Gearbox ratio 


C, : Power coefficient 

P „ : The aerodynamic power 
T,, : Torque on the high speed 
T, : Torque on the low speed 
R: Radius of the wind turbine rotor 
v, : Wind speed 

v,: The blade tip speed 

A : Tip speed ratio 

p : Air density 

r, : Stator resistance 

x, : Stator leakage reactance 

: Slip 

x,, : Magnetization reactance 


nA 


q 
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Xp, : Rotor leakage reactance 
rp, : Rotor reactance 
I. INTRODUCTION 


Wind energy is an energy source used for centuries old: it is 
a renewable energy source that saves fuel which contributes to 
the protection of the environment, non-polluting and 
economical. It is a source of inspiration for many 
manufacturers; they always push them to think new robust 
solutions, allowing maximum energy extraction and study of 
systems feasible. 


Currently, Tunisia exploits its potential energy, including 
wind energy. The power and the availability of the wind are 
random and depend on the site and the period considered and 
Tunisia considers moderately windy [8], [4]. 


The STEG (Tunisian Society of Electricity and Gas) 
integrated the wind one among its choices in the objective to 
contribute to preserve the environment. It engaged of the 
feasibility studies for the exploitation of wind engage in the 
production of electricity, in cooperation with company MADE, 
carried out a first wind power station with Sidi-Daoud whose 
commissioning of the first slice took place in August 2000, [7]. 


In this paper, we present in the first section, the wind power 
farms of Sidi-Daoud and we describes the different models of 
the wind turbines constitute this park. In second section, we 
will perform simulations of part of this Wind farms. First of all, 
for the validation of the results, a software PSAT dedicated to 
the study of the Power Systems was introduced especially that 
this software is interested in the study of the wind systems. The 
last section of this paper present the simulation results of the 
wind farm under the turbulent wind upstream of each rotor 
turbine. 


Il. WIND ROSE OF THE WIND FARM OF SIDI- 
DAOUD 


The wind rose is a spatial representation of the variation of 
the wind direction and average speed for such a Park. 
It illustrates the direction of the dominant winds ona site and 
allows planning the installation of turbines in order 
to minimize the effect of wake caused by the neighboring 
obstacles. 


Fig.1 represents the wind roses of wind farm of Sidi-Daoud 
with 36 directions. We note the importance of the wind coming 
from the sectors West and South-east, and the existence non- 
negligible wind from the sector northwest. Furthermore, the 
wind was calm persists both sectors Southeast and South [2]. 
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Fig. 1.Wind Rose in the region of Sidi-Daoud [6] 
Figure 2 presents the frequency of wind speeds at the park 


since the 01/01/2010 until 31/12/2010 the measured values 
represents the average value of the wind speed for hour. 
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Fig 2. Frequency of wind speeds from 1/01/2010 to 31/12/2010 in the wind 
farm of SidiDaoud [14]. 
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UI. PRESENTATION OF THE WIND POWER FARM OF SIDI- 
DAOUD 


The total surface used for the establishment of the power 
plant is about 9 ha and extends on 3.5 km from the coast shown 
in Fig.3 It is about a site sufficiently been windy, capable to 
receive and close to the grid, [5]. 


Wind turbines of the first and second section are 
established on the summits of the two mountains "Djebel El 
Hammam" and "Jebel Ghormane" whose altitude is 
respectively 50 and 100 m above the sea level. The wind 
turbines of the third section are located at the bottom of these 
two hills and approximately one hundred meters of marine 


coasts [2]. Figure 3 shows the topographic profile and the 
location of wind turbines on site. 
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Fig. 3.SidiDaoud site relief [2] 


IV. DIFFERENT MODELS OF WIND TURBINES CONSTITUTE 
THE CENTRAL OF SIDI-DAOUD 


The wind farm of Sidi-Daoud consists of four models of 
aero-generator with the horizontal axis, all from the same 
manufacturer Spanish MADE. The first model was installed is 
the AE32 in 2000 with a power of 330kW, using an 
asynchronous squirrel-cage generator and a system of 
aerodynamic unhooking “stall” with a fixed pitch angle of the 
blades. The site comprises 32 aero-generators of this model 
“phase A”. 


A first extension of the site in 2004 is then accomplished 
with the implantation of 12 wind turbines "Phase B" with 3 
different models, 10 AE46 whose power of each one is of 660 
kW using same technology as its predecessor, the AE52 with a 
power of 800kW using a synchronous generator, with a 
variable pitch angle of blade and finally the AE61 whose 
power is of 1320 kW and using the same technology as the 
AE32 and the AE46. 


Following the good wind potential of the site and to the 
infrastructure already exists, a second extension "phase C" 
was to realize in 2009 comprising 26 AE61 generator, we 


summarizes the three stages of wind farm of Sidi-Doud in a 
TABLE I. 


TABLE I. TABLESTEPS SIDI DAOUD WIND FARM [8], [2] 
Step A first second 
extension | extension 
Number of turbines 
installed power (Mw) 10.56 34.32 


designations 


Industrial August | September 
constructor MADA MADA MADA 
| constructor T | Spain 
(GWh) 





The STEG is provided the characteristic of the power at the 
level of the turbine according of the wind speed for each aero- 
generator of the central Sidi-Daoud. This characteristic 
illustrated by Fig.4. The machines start from the same speed of 
3 m/s (except AE-32, which starts at 4 m / s) this speed is boot 


speed from which the turbines begin to produce and must stop 
at maximum speed 25 m/s. 


This power reaches its nominal value for a nominal wind 
speed of about (13, 15, 12 and 17 m/s). Beyond the nominal 
speed, the power supplied by the synchronous machine AE-52 
remains constant, on the other hand, that provided by 
asynchronous machines AE-32, AE-46 and AE-61 decreases 
slightly with the wind speed [2]. 
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Fig. 4. Power curves of wind turbines 


V. WIND TURBINE MODELING 
A. Aerodynamic model: 

The mechanical P, power available on the shaft of a wind 
turbine v, of a captured any wind speed is determined 
analytically by the following formula: 

pe + oC, (A) AV (1) 
2 


W 


With A. is the area swept by the wind turbine (m°) and 


p is the density of air (kgm? ), whose value depends on the 
height at which is installed the turbine and the power 
coefficient C, depends on the speed ratio which is defined as 


follows: 


1 
100 = 


L 


With: 4 =— — (3) 


| 0.002 
A 


The speed tip ratio 2 is the ratio between the blade tip 
speed v, and the wind upstream the rotor v, : 


v 2Rw,,, 
A=— =Nep 


Vy PVw 





(4) 


B. Transmission System Model 


We can model the driving device of the chain conversion of 
wind turbine by a mechanical model approximate with two 
mass shown in fig.5[3],[1]. 
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Fig. 5. Schema of the transmission system with two masses [3] 


The differential equations of the mechanical part by a 
model of a two-mass drive train model obtained by neglecting 
the coefficients of external frictions of the rotor and generator 
are given by: 


AW yy _ (Tor —K, y) 


dt 2H 
dw, (K,y-T,) 
Se a (4) 
dt 2H, 
y= Q; (Wyr = Wm) 
With T, „ is the mechanical torque defined by: 
F 
T, = — (5) 
wW 


C. Induction Generator Model 


The simplified electrical circuit used for the induction 
generator with squirrel cage is the same as that for the 
induction motor with single cage, shown in Fig.6. The only 
difference is that the currents are positive when injected into 
the network [10]. 


r x XI 





Fig. 6. Order III induction motor: electrical circuit [10]. 


The equations are formulated in terms of the real (r) and 
imaginary (m) axis, with respect to the network reference 
angle. In a synchronously rotating reference frame, the link 
between the network and the stator machine voltages is as 
follows: 


v, =V sin(—0) 


(6) 
Vm =V cos(@) 
And the power absorptions are: 
P=v i +V by 
(7) 


men od . P29 
Q=V,,1, —Vrly +b, (v; +va ) 


Where be is the fixed capacitor conductance which is 
determined at the initialization step. The differential equations 
in terms of the voltage behind the stator. 


The link between voltages, currents and state variables is as 
follows: 


e, —Vv, =ri,—xi 


i (8) 
CVT 
and: v, =-r,i,„ + xi, ag e, (9) 


Where e, is the real part of transient voltage, e, is the 


transient voltage on imaginary axis and x is transient 
reactance. 


(10) 


Where x), x and T can be obtained from the generator 
parameters: 


S m 
Xp X 
rea (11) 
XR] FX 
' Xp F Xn 
‘0 Q 
bĪRI 


Where the electrical torque T, is: 
T, =e i, +e i, (12) 
VI. PV GENERATOR 


PV generators fix the voltage magnitude and the power 
injected at the buses where they are connected, as follows [10]: 


P=-P, 
(13) 


VII. WIND MODELS 


On PSAT the wind is modeled by a block window allows to 
easily analyze and modify all the parameters related to the 
behavior of the wind as the wind model type as the Weibull 
distribution or the composite model which includes average 
speed , the intensity of wind gusts and turbulence level or by 
real measures of the wind speed [3]. 


In our case we chose the method most frequently used to 
describe the wind speed which is the Weibull distribution. 


The variations in wind speed a function of time v, (t) are 
obtained by distributing wind by the equation (3.1). 


Vy (t) = Gel (14) 


C 


Where q (t)a generator of random numbers, C is the scale 


factor of the Weibull distribution and K is the shape factor. 
Then, the wind speed is calculated, in framing the initial 
average value of the wind speed V,, as follows [10]: 


v(t) = (+v, (t)— Cn Vina (15) 
Where Ŷ, is the mean value ofv,, (t). 


VII. SIMULATION OF PHASE C OF THE SIDI DAOUD WIND 
FRAM 


The phase C of Sidi-Daoud wind farm consists of 26 aero- 
generators AE61 spread over 3 lines, the first line designated" 
1C" contains 9 aero-generators located about 1.56 km from the 
substation, the second line contains 10 aero-generators located 
about 2.765 km from the post of transformation and the third 
line is located at 0.77 km from the post of transformation and 
contains 7 aero-generators. In Sidi-Daoud wind farm, 
converting electricity produced in two stages, the first of 690 V 
to 30 kV by a transformer of the each turbine, and the second 
from 30 kV to 90 kV by the external post after passage of the 
MT voltage underground cable from the foot of the tower to 
the MT / HT [11]. We visualize the wind turbines in real 
operating and identify the main features. The simulated 
structure is shown in fig.7.Parameters of the aero-generator 
constituents AE-61 are recapitulated in Table II. 
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Fig. 7. Phase C of the Wind Farm of Sidi Daoued 
TABLEII.PARAMETERS OF THE WIND TURBINE COMPONENTS TABLEIII. THE RESULTS OF ACTIVE AND REACTIVE POWER IN THE 


PERMANENT REGIME FOR PHASE C 


a ee lee ee es) (pu) (pu) 


ee a 
[3 [ es p 0.1858 | one | 
o4 | 85 | 03587 | -0.2329 | 
PST S| o O 0.1519 
6 | 57 | 01032 | 0.06742 | 
7 | 66 | 01724 | -0.01706 | 
8 | 64 | 01557 | 0.00275 | 
po ft 43002995 | 0.15743 

PSAT we provided in the first time, after resolution of the 
equations powers a static report that contains the results of the 
various state variables of the system, this part is a step to 
initialize the wind turbines. Generators PV are necessary to 
impose the desired voltage and active power on the bus of wind p24 | 6 | 01245 | 004256 
turbines. We extracted the static results of active and reactive 


power in the bus of each machine. These results are 
summarized in the following Table ITI: 


—— values 


| 
Ll 
— 
a 





Follows, PSAT can also perform the simulation in the time 
domain in order to study the behavior of the dynamic system. 
These simulations are based on numerical integration of 
differential algebraic equations (DAE) that describe the 
dynamic function of the system in our case they are used to 
study the effect of variation of wind speed on all variables. 


The 26 wind turbines of Sidi-Daoud farm are simulated for 
turbulent wind speed with a scale factor of the Weibull 
distribution c = 2 and a shape factor k = 9. Figure 8 shows the 
variation of the wind speed for wind turbine number 1. We 
note that the wind speed oscillates around its average value 7.1 
m/s. 
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Fig. 8. Wind speed 


Figures from 8 to 15 shows that the behaviors of all 
variables are random and oscillates around its operating values, 
it helps to explain that these variables vary due to the variation 
of wind speed. 
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Fig. 9. Active power in the bus network 
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Fig. 10. Reactive power supplied by the network 


Figure 8 and Fig.9 shows respectively the active power and 
reactive power in the bus which is connected to the network. 
Gaits powers are random and fluctuate around its average 
values. We note that the network received total power of 4.2 pu 
is equivalent to a production of 6.5 MW and total reactive 
power supplied around 1.9 pu which is also equal to 2.9 
MVAR, which gives a power factor equal to 0.91. 
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Fig. 10. Active power in the 3 lines 1C, 2C and 3C 
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Fig. 11. Reactive power in 3 line 1C, 2C and 3C 


Figure 10 and Fig.11 illustrate respectively the active and 
reactive power in different lines 1C, 2C and 3C. Three lines 
produce successively 1.7, 1.17and 1.25pu active powers of the 
network and received respectively 0.29, 0.45 and 0.13 of 
reactive power. 
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Fig. 12. Real component of the fem for two different values of wind speed 
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Fig. 13. Imaginaire component of the fem for two different values of wind 
speed 
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Fig. 14: currents for two different values of wind speed 
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Fig. 15. power of the machine for different values of the two wind speed 


We shows the real and imaginary component of the fem of 
the machine for two values of wind speed 8.5 m / s and 4.3 m / 
s Fig.12 and Fig.13. We notice that the magnitude of the 
transient f.e.m increases when wind speed increases, also for 
currents and voltages in the bus machines Fig. 14 and Fig. 15. 
We conclude that the variation of the wind speed value affected 
the production of Wind Energy. 


VX. CONCLUSION 


This paper presents the simulation of phase C of Sidi- 
Daoud wind farm composed by 26 wind turbines at constant 


speed using such AE61gerator directly connected to the electric 
grid. The simulation aims to better understand the dynamic 
behavior and the electrical properties of the constant speed 
wind turbine. In this work, we demonstrate that the electrical 
grid is affected by the turbulent wind speed variations it can 
degrade the quality of power and has a detrimental effect on 
the installation. The greatest handicap of the constant speed 
wind turbine technology studied (AE61) is that it does not 
offer the possibility of the power control under wind speed 
variation. 
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